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Liquid markets and market liquids 

Collective and single-asset dynamics in financial markets 

Gianaurelio Cuniberti a and Lorenzo Matassini 

Max-Planck-Institut fur Physik komplexer Systeme, Nothnitzer Strafie 38, D-Of 187 Dresden 
Received: / Revised version: 

Abstract. We characterize the collective phenomena of a liquid market. By interpreting the behavior of 
a no-arbitrage N asset market in terms of a particle system scenario, (thermo) dynamical-like properties 
can be extracted from the asset kinetics. In this scheme the mechanisms of the particle interaction can be 
widely investigated. We test the verisimilitude of our construction on two-decade stock market daily data 
(DAX30) and show the result obtained for the interaction potential among asset pairs. 

PACS. 02.50.Sk Multivariate analysis - 89.90.+n Other areas of general interest to physicists 



Since the late 80s, with the introduction of electronic trad- 
ing, huge quantities of financial data became available (or 
at least on sale) for both investment and research analy- 
sis. Quite unusually outside the natural science panorama, 
this novelty opened the way to test the reliability of theo- 
ries and conjectures about the behavior of financial mar- 
kets. One of them, a paradigm for financial mathemat- 
ics, is the random character of markets [l|, that is un- 
predictability. It has recently been proved, nevertheless, 
that a certain degree of correlation is still present on ex- 
tremely short time scales Q]. Despite that, the intermedi- 
ate scales are dominated by random behavior with Levy 
stable statistics of asset returns Q . The possibility to ex- 
tract information on the future evolution of a single asset 
by knowing a big enough ensemble of its past values mat- 
ters indeed institutional traders, who can generally inter- 
vene on the market in real time (with delays smaller than 
few seconds). Their presence reduces at minimum time 
correlations and consequently speculation possibilities. 
Time dependence is however only one possible domain 
for surveying similar patterns inside financial signals. The 
other domain for correlation detection, whose exploration 
was greatly facilitated by modern computation facilities 
is the 'spatial' one. In fact, albeit much efforts are spent 
in studyingcorrelations in the time dynamics of a single 
asset (see Q and || for a digest of the recent physicist 
and economist approach, respectively), there are many 
applicative and fundamental reasons for understanding 
deeply spatial, commonly referred as multivariate, corre- 
lations. A financial market is not simply a juxtaposition 
of different prices which are organised on an independent 
basis, but rather a complex system of interacting con- 
stituents 0. The latter are then monitored by sampling 
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single prices with respect to an arbitrary currency. Hence 
the study of correlations among different asset time sig- 
nals is of peculiar importance. By the way, this is also 
the case in many problems involved in the modern risk 
management theory, where the composition of a certain 
portfolio strongly depends on the movements of different 
underlying assets. On a more fundamental level, the in- 
teresting issue is the comprehension of how price changes 
can be separated, with a sufficient degree of confidence, in 
single asset- and collective- behavior. 
Since the Markowitz's work on the theory of optimal port- 
folio 0, much effort has been spent to characterize cor- 
relation matrices of financial assets JsJ . In recent contri- 
butions, different physics concepts have been adopted to 
endeavor this type of problem, mainly because the study 
of correlations represents a paradigm of a wide class of 
physical problems for which powerful tools have been de- 
veloped. A bivariate analysis of the futures on the German 
and Italian bonds showed that despite the perfect uncor- 
relation of the single tracks, the crosscorrelation of the two 
signals was significantly non zero: the signals considered 
described two random, but similar, processes This be- 
havior emerges quite generally in the stock market, where 
certain asset clusters 'move' in a particularly correlated 
way with respect to remaining titles. Using equal time 
cross-correlation matrices and several physics-borrowed 
tools such as the random matrix theory, these conjectures 
have been quantified p0| . In a recent study, the structure 
of a N stock market has been investigated as regarding 
the multivariate structure in a global window period [fj"l| . 
In this paper, we propose a method to investigate asset 
correlations by interpreting asset growth rates as observ- 
ables of a particle system scenario. This idea is carried out 
by introducing a formal map between the logarithmic re- 
turns and the distances among gas particles. The strength 
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of this analogy resides in the possibility to separate col- 
lective motion from the single asset dynamics through the 
investigation of mutual interactions among titles. Wielded 
by the theory of liquids, we can study the thermodynam- 
ics of the system and interprete its temperature as a mea- 
sure of spatial volatility, as compared with the more famil- 
iar (temporal) volatility. The 2-asset interacting potential 
is then calculated on the isothermal (isovolatilc) market. 
In the remainder of this paper a time dependent asset- 
distance and a moving frame model are introduces. The 
implementation of this scheme is performed on daily stock 
market data taken among the 30 most capitalized titles 
forming the Deutscher Aktien indeX (DAX30) in the pe- 
riod 30 Dec 1987 to 7 Mar 1995 (1800 trading days). To 
maintain a continuity of quotation, we have selected the 
maximal subset of 23 assets which, in the above mentioned 
period, remained in the DAX30 basket and did not split. 
Our discussion and comments conclude the paper. 
As a general starting point, we consider a collection of 
asset, which is a suitable subpart of titles in a stock market 
(better if one representative for every economic sector), a 
collection of currency prices, or any combination of them. 
The value of the asset fli at time t, is expressed in unity 
of asset fij by means of conversion factors Pij (t) : 



fli(t) 



(i) 



The indices i and j span all N considered assets form- 
ing the market. By writing eq. (0) for another couple of 
indices, a no-arbitrage equation for a liquid market is ob- 
tained Pij = PikPkj- Its multiplicative symmetry is re- 
flected in a corresponding additive symmetry of the loga- 
rithmic returns 



dUt) = - log 



Pij(t-r a ) 



(2) 



where r a <H is a collection of H time horizons. The rescal- 
ing of the log-returns to the considered time horizon is 
solicited by its interpretation; in the idealized limit of 
prices with (deterministic) growth laws, we get Pij(t) oc 
exp (dijt), so that the quantity defined in eq. (0) turns out 
to be the growth rate between asset i and j, independently 
on the time horizon. The latter can be considered as a 
long term limit when one refers — for example — to prices 
of stocks with respect to currencies. In the opposite limit 
of extremely small returns (which eventually corresponds 
to short time lags), d a is the rate of the absolute return, 
d = AP/(PAt), obtained by logarithmic expansion. 
As position ^) points out, the display of the time series Pij 
by arranging them in the H dimensional variable , gives 
a natural embedding for a dynamical system oriented anal- 
ysis 12 . This is not difficult to understand when thinking 
that the log-return on a certain time horizon r* is propor- 
tional to the average of log-returns on sub-multiples of t* . 
Thus the component d a can be written as a linear com- 
bination of delayed components d a<a . The no-arbitrage 
symmetry of the log-returns hints for the further identi 
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Fig. 1. Time dependence of the correlated volatility a, and 
the temperatures T (shifted as a visual aid) relative to the x— 
and r-coordinates (a); the corresponding PDFs are in panel b. 
All the calculations refer to four horizons (H — 4) of 1, 5, 20, 
and 250 market days. 



mcnsional euclidean space (in the present work we choose 
the canonical one) induces a well defined distance^] ||dy|| 
between asset i and asset j. As an intrinsic character of 
financial markets no asset can be regarded a priori as an 
absolute quantity, that is why we ended up only with mu- 
tual distances among asset. Nevertheless some truly single 
asset property can be extracted by the symmetry of the 
problem and interpreted consequently The matrix is 
skew symmetric, ergo diagonalizablc; its spectrum is en- 
tirely on the imaginary axe JL3[ . One of its three different 
eigenvalues is zero and the corresponding eigenspace is or- 
thogonal to 1 =(1,1,..., 1)* and it = (xi, X2, . . . , x^)', 
where 



1 N 



(•3) 



j=i 



here the arrows indicate super-vectors (vectors in a H x N 
space). The two remaining eigenvalues are ±iNa corre- 
sponding to the eigenvectors 1 =F vSt /a, where 



E 



(4) 



l<i<j<N 



n ,. r j / j\ j- , i In particular let us observe that x, —x, = d,,-., , to say that 
hcation of dy as an (oriented) distance vector between _ 1 3 y ' J 



asset i and j; in fact (a.) An = 0, (6.) dy = — dji, (c.) 
dij = dik + dfcj . It is easy to see that any norm in a H di- 



1 Here the three distance-defining axioms are obtained sim- 
ply by properties a.-c. 
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Fig. 2. Masses rm as calculated after eq. (g) versus the as- 
set label |14| . Crosses indicate a;-frame calculations with mean 
1.195, circles the r-frames with mean 1.184. 



we have introduced a frame in which every single asset 
is assigned to an absolute position: the problem of the 
behavior of the N assets of the market is now translated 
to a physical problem of N interacting particles (a liquid) 
in H dimensions, with coordinates Xi, X2, . . . , Xjv. At time 
t, Xi(t) is the H dimensional position of particle i. Note 
that, according to its definition, the distance between two 
assets is zero when the price of one with respect to the 
other remains constant. Furthermore, it is easy to check 
that the x's vectors are centered, hence the positions x, 
are referred to a coordinate frame which attributes to the 
center of mass of our liquid a trivial dynamics. From the 
financial point of view, it states the closure of our system: 
the N assets are watched as complementary, with zero 
overall return. This does not mean that the applicability 
of the present construction is restricted to those market 
where this property is nearly fulfilled (as an example in 
the foreign exchange). In stock markets, which experience 
escape and retention events that is positive and negative 
return periods, the x are automatically selected within a 
neutral frame which keeps track of the particle cloud. Of 
course nothing prevents from starting the analysis of an 
extended market with a huge number of constituent assets. 
Some of them would follow similar dynamics by evolving 
in a closer cluster with respect to others. This could help 
in order to reduce N to a lower number without losing the 



basic features of the liquid behavior |15 
Coming back to the map construction it is easy to show 
that as a consequence of the centered character of the x's 
coordinates, a is exactly their standard deviation. Its Hth 
power is a measure of the volume of our system. The finan- 
cial counterpart of it is what we call correlated volatility, 
so to stress that it is a quantity merely connected to the 
spatial interactions of the particles at a certain time. As 
the usual volatility takes into account the temporal vari- 
ability of an analyzed fixed asset, we are here referring to 
a measure of a spatial variability of a group of interacting 




Fig. 3. (a) Plot of the pair potential u(r) for the whole data 
set over four horizons (H — 4) of 1, 5, 20, and 250 market 
days, and (b) the time distribution of the inter-asset distances 
(a better resolution figure is available upon request). 



assets at a fixed time. Moreover, even after the compen- 
sation of the split discontinuities, the correlated volat ility 
shows clusterization around bubble and crash periods []15| . 
We take now advantage of the structure of the eigenvec- 
tors of the distance matrix and rescale the x coordinates to 
volume renormalizcd ones = Xj/<7, their difference is ac- 
cordingly Yi—Yj = dij/a, so that the two non trivial eigen- 
vectors are indeed 1 =Fi~r\ The r-frame, being the solution 
of the eigenvalue problem for the distance matrix, is a vol- 
ume preserving frame. Once the volume of the system is 
stabilized, one may wonder which is the dependence of the 
liquid temperature on time. Thus, by analyzing the empir- 
ical behavior of the ensambled averaged square (finite dif- 
ference) velocities v, (i) = (r$(t) — rj(t — Ti ))/ti , we found 
that the r-system is thcrmostatcd at a fix temperature 
T = ((vf(t)^ \ /H ; the correlated volatility is therefore 
a measure of the temperature of our system. Fig. [l] shows 
this fact: in panel (b), it is plotted the time dependence 
of the correlated volatility a and of the temperatures T[x] 
and T\r] calculated by averaging the square velocities in 
the x— and r-frame, respectively. In order to contrast the 
results, the time averages of a and T[x] are rescaled to 
T (the time average of T[r}). The scale of T is in fact 
fixed by the underlying assumption of an unitary Boltz- 
mann's constant. To check possible ergodicity properties 
of the system we have also analyzed the time averaged 
square velocities of the single assets and extracted from 
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them mass terms 



HT 



(5) 



Fig. H shows that the masses are only slightly affected by 
the reference frame used to calculate them. This indicates 
that they are an intrinsic property of the asset regardless 
of the kinetics details. To prove this statement, we have 
plotted, in Fig. the correction to the temperature due to 
the asset masses T m [r] = (m,ivf(t)') i / H. In order to inves- 
tigate the nature of the interaction of the particle system 
under study, we have calculated the two point correlation 
function lfl6| 



9(r) 



N(N 



YK § (r-\Mt)-r j (t)\\)) t 



and the related pair potential u(r) oc — log g(r). In Fig. ^|., 
the potential u(r) is shown. The great distance tail of u(r) 
is linear (correlation coefficient= 0.9994, for a regression 
in the region 2 < r < 4 over 446 points giving the line 
u = ar + b, with a = 0.689±0.001 and b = -1.101 ±0.004) 
indicating the strong long range attraction of the market 
liquid. On the other hand at small distances two different 
behaviors emerge. By decreasing the asset-asset distance 
an equilibrium point is reached. At smaller distances a 
barrier is presented, followed by a region corresponding to 
less intense repulsive forces. We interprete it as a signature 
of the inhomogeneity of the system, which allow at small 
distances the formation of privileged pairs (clusters). As 
a consequence, we expect that in a wider market (here we 
consider the quite diversified but small pool of the DAX30 
assets) this tendency could even be more pronounced. 
To conclude, we have introduced an interpretation scheme 
for the returns of a N asset market. Here, the validity of a 
no-arbitrage condition is guaranteed by the assumed liq- 
uid character of the market. By implementing the embed- 
ding, naturally prompted by the structure of the returns, 
we have been able to map the financial signals in positions 
of particles of an interacting gas (a liquid). Therefore by 
the only means of the geometrical construction, we have 
given purport to the eigenvalues and eigenvectors of the 
distance matrix diagonalization problem as the correlated 
volatility and thermostated coordinates respectively. One 
of the strength points of this method is its easy general- 
izability to the case of great TV, albeit here we have re- 
stricted our analysis to a relative small asset market. 
On the other hand a word of caution is needed in a great 
N market. The results presented here share the plain as- 
sumptions of isotropy and homogeneity of the market liq- 
uid. Indeed they should become weaker for very large 
and differentiated markets. There, the pair potential in- 
troduced here is supposed to maintain the same great- 
distance properties (linearity). At the low-distances (where 
clustering emerges), in analogy to what is done in the 
study of ionic liquids Jl7| ], a generalized pair potential 
could be introduced in order to include both anisotropy, 
cluster formation, and specie diversification: these issues 
are under investigation and will be published elsewhere fig] . 



Besides, this approach is straightforwardly employable for 
time dependent clustering. A procedure similar to the one 
adopted to organize static distances between assets in hi- 
erarchy trees, given in Ref. [O], could be generalized to 
the time dependent distance matrix (^). 
From the financial perspective, the construction presented 
in this paper is following a sort of TEsop's the fox and 
the grapes strategy. It is easy to despise what one cannot 
get and in quantitative finance the scarce goods are the 
rare events. Since there is no methodology to deal with 
misprediction given by the insufficiency of statistics, we 
try to wash it out from the dynamics by exploiting the 
symmetries of the problem. After all, as we have shown, 
the calm (non-bubble, non-crash) side of financial markets 
has anyhow a lot to say. 

The data for the empirical analysis were kindly provided 
by Deutsche Borse AG. We would like to acknowledge 
fruitful discussions with A. Amici, F. Lillo, R. Mantegna, 
E. Scalas, and U. Tartaglino. 
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